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High-enegy effective action from scattering of QCD shock waves 
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At high energies, the relevant degrees of freedom are Wilson lines - infinite gauge links ordered 
along straight lines collinear to the velocities of colliding particles. The effective action for these 
Wilson lines is determined by the scattering of QCD shock waves. I develop the symmetric expansion 
of the effective action in powers of strength of one of the shock waves and calculate the leading term 
of the series. The corresponding first-order effective action, symmetric with respect to projectile 
and target, includes both up and down fan diagrams and pomeron loops. 
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I. INTRODUCTION 

It is widely believed that the relevant degrees of free- 
dom for the description of high-energy scattering in QCD 
are Wilson lines - infinite straight-line gauge factors. An 
argument in favor of this goes as follows pj . As a reslut 
of a high-energy collision, we have a shower of produced 
particles in the whole range of rapidity between the tar- 
get and the spectator. Let us demonstrate that the in- 
teraction of gluons with a different rapidity is described 
in terms of Wilson lines. Consider the fast particle inter- 
acting with some slow gluons. This particle moves along 
its classical trajectory - a straight line collinear to the ve- 
locity, and the only effect of the slow gluons is the phase 
factor Pexp{ig J dx^A^} ordered along the straight-line 
classical path (here describes the slow gluons). This 
picture is reciprocal - in the rest frame of fast particles 
the fast and slow gluons trade places: former slow gluons 
move very fast so their propagator reduces to a Wilson 
line made from the (former) fast gluons. We see that the 
particles with different rapidities perceive each other as 
Wilson lines and therefore these lines must be the rele- 
vant degrees of freedom for high-energy scattering. The 
goal of this approach is to rewrite the original functional 
integral over gluons (and quarks) as a 2 + 1 theory with 
the effective action written in terms of the dynamical 
Wilson lines. 

For a given interval of rapidity, the effective action is an 
amplitude of scattering of two QCD shock waves, see Fig. 
n Indeed, let us integrate over the gluons in this interval 
of rapidity 771 > 77 > 772 leaving the gluons with r\ > r\\ 
(the "right-movers') and with r\ < 772 (the "left-movers') 
intact (to be integrated over later). Due to the Lorentz 
contraction, the left-moving and the right-moving gluons 
shrink to the two gluon "pancakes" or shock waves. The 
result of the integration over the rapidities 771 > 77 > 772 
is the effective action which depends on the Wilson lines 
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FIG. 1: High-energy effective action as an amplitude of the 
collision of two shock waves. 



made from the left-and right-movers. 

Due the parton saturation at high energies the 
characteristic scale of the transverse momenta in hadron- 
hadron collisions is Q s ~ e CI? [U EJ Q and therefore 
the collision of QCD shock waves can be treated using 
semiclassical methods 0. Within the semiclassical ap- 
proach, the problem of scatering of two shock waves can 
be reduced to the solution of classical YM equations with 
sources being the shock waves E3 (see also [llj). At 
present, these equations have not been solved. There are 
two approaches discussed in current literature: numeri- 
cal simulations [12J and expansion in the strength of one 
of the shock waves. The collision of a weak and a strong 
shock waves corresponds to the deep inelastic scattering 
from a nucleus (and scattering of two strong shock waves 
describes a nucleus- nucleus collision). The first term of 
the expansion of the strength of one of the waves was cal- 
culated in a number of papers 

111 Ellli- Recently, the 
classical field was calculated up to the second order in a 
weak source 0] . I will use some formulas of Ref . 0] , al- 
though the main result for the effective action will be de- 
rived independently. The obtained effective action, sym- 
metric with respect to projectile and target, has a "built- 
in" projectile-target duality (which is a highly nontrivial 
property of the light-cone Hamiltonian in the framework 
of the Hamiltonian approach |l7lllclll9ll2Ctl2l| ). In terms 
of Feynman diagrams the effective action includes both 
"up" and "down" fan diagrams and therefore it describes 
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FIG. 2: Rapidity factorization. 



pomeron loops which are a topic of intensive discussion 
in the current literature fTl HI HH I2H l2lll2^| . 

The paper is organized as follows. Sec. 2 is devoted to 
the rapidity factorization which is the starting point of 
the shock- wave approach. In Sect. 3 I define the high- 
energy effective action as a scattering amplitude of QCD 
shock waves and develop the expansion in commutators 
of Wilson lines. In Sec. 4 I find the effective action for a 
given (infinitesimal) range of rapidity in the leading order 
in this expansion. The corresponding functional integral 
over the dynamical Wilson-line variables is constructed in 
Sec. 5. The explicit form of the first-order classical fields 
created by the collision of two shock waves is presented 
in the Appendix. 



II. RAPIDITY FACTORIZATION 

The main technical tool of the shock-wave approach 
to the high-e nergy sc attering is the rapidity factorization 
developed in jllL l23j . Consider a functional integral for 
the typical scattering amplitude 



DAJ(p A )J(p B )J(-p' A )J(-p' B ) 



(1) 



where the currents J (pa) and J(pb) describe the two 
colliding particles (say, photons). 

Throughout the paper, we use Sudakov variables 



k = apt + /3p 2 + fcj 



and the notations 



x 9 — Py x^ 



P2 x ^1 



2 X ' X 

s - 

—x, x 
2 ' 



1 
1 

71 



(2) 

(x° - x 3 ) 

(x° + x 3 ) (3) 



Here p\ and p 2 are the light-like vectors close to pa and 

Pb- Pa=Pi + E fP2, Pb = Pi + Pi- 
Let us take some "rapidity divide" 771 such that t\a > 
Vi > Vb and integrate first over the gluons with the ra- 
pidity 77 > 771 , see Fig. |2 From the viewpont of such 
particles, the fields with 77 < 771 shrink to a shock wave 
so the result of the integration is presented by Feynman 





FIG. 3: Shock wave in the temporal gauge. 



diagrams in the shock-wave background. With the LLA 
accuracy, in the Feynman integrals over the gluons with 
77 > 771 one can set 771 — > — 00 (replace the "rapidity 



divide" vector ei = pi 



Pi by the light-like vec- 



tor pi) so the shock wave is infinitely thin and light- 
like. In the covariant gauge, this shock-wave has the only 
non-vanishing component A. which is concentrated near 
x* = 0. In order to write down factorization we need to 
rewrite the shock wave in the temporal gauge A = 0. In 
such gauge the most general form of a shock-wave back- 
ground is (see Fig. |3J) 



A i =U{6(x*) +l40{-x*), A. = A* = (4) 



where 



Ut-diUi, 

g 



w, 



U\-diU 2 
9 



(5) 



are the pure gauge fields (filling the half-spaces x* < 
and x* < ). There is a redundant gauge symmetry 

Uxixx) -> Ux(x x )to{xx), U 2 (xx) -» tMzxMzx) 

(6) 

related to the fact that gauge invariant objects like the 
color dipole 



Tr{[oop 2 , -OOp2]:r 

x [y ± + oop 2 ,x_i_ 



[xx -oop 2 ,y± - oop 2 }[-oop2,oop 2 }y 
f oo77 2 ]} ~ Tr{U lx UiU 2y ul y } (7) 



depend only on the product U\ z U 2z . In papers 0, |n| 
this symmetry was used to gauge away U 2 and simplify 
the shock wave to Ai = Ui6(x*) while in Ref. [l6( the 
opposite case U± = (Ai — Ui6(—x*)) was considered. 
In the present paper we keep this gauge freedom - as we 
shall see below it simplifies the effective action for the 
Wilson-line integral. 

The generating functional for the Green functions in 
the Eq. J3} has the form (cf. 0) 



DAJ(pa)J(-p' a ) e lS ^ + ^ d2z ^ F " fi) ^ u "~ u "^ (8) 



where (F el = e^F^ etc.) 

I' OC 

(0,F ei ,0) z = du [0,ue] z F ei (ue + z±)[ue,0] 2 



d 

[O,ooe] x (i^-T + gAi(ooe + z±))[ooe,0] z (9) 
d 

[0, -ooe] z (i— t ■ + gAi(-ooe + z±))[-ooe, 0] z 
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and (0,F Ml ,0) a = 2tr t a (0, F^, 0). (Throughout the pa- 
per, the sum over the Latin indices runs over the 
two transverse components while the sum over Greek in- 
dices runs over the four components as usual). 

It is easy to see that the functional integral 10 gener- 
ates Green functions in the Eq. (@J background. Indeed, 
let us choose the gauge A* = for simplicity. In this 
gauge, (0, F* i: 0) a = Ai(oop 2 + z±) - A t (-oop 2 + z±) so 
the functional integral JSJ takes the form 



DA J( PA )J(-p' A ) 



(10) 



x e 



iS(A)+i[ d 2 z ± {A i {oop2+z ± )-A i (-oop 2 +z±)) a {U^ i -U^ i ) z 



Let us now shift the fields A{ — > A{ + Ai and where 
A 1 = a;*). The only non-zero components 

of the classical field strength in our case are F,i = (Uu — 
14h)5(jX*) so we get 



S(A + A) = - [ d i zD i F i .A* - - I d 2 z ± dz. 

8 J S 



x A l F mi 



±A»(D 2 9fl „ - 2%F^)A V + ... (11) 



In the A* = gauge the first term in the r.h.s. of Eq. 
(|1 1|1 vanishes while the second term cancels with the cor- 
responding contribution ~ — (Ai(oop2 + z±) — Ai{— oop2 + 
z±)) a U ai coming from the source in Eq. (JSJ. We obtain 



DA J( PA )J(-p' A ) (12) 

x e iS(.A)+iJ d 2 z ± (A,(oc P 2+z ± )-A t (-cop2+z ± )) a (U? i -U^)z 

= J DA J(j>a)J(- ■ 



n' \ASd 2 zA^{D 2 gfi y-2tF^)A" 



which gives the Green functions in the Eq. back- 
ground. 

To complete the factorization formula one needs to in- 
tegrate over the remaining fields with rapidities r) < rji: 

J DADAJ( PA )J( PB )e- lS ^ e iS ^J[-p A )J{- PB ) 
= J DAJ(p A )J{-p' A ) J DBJ{p B )J(-p' B ) 



x e 



iS(A)+iS(B)+if d 2 z x (0,F eii ,0)J(0,G ei i,0)° 



(13) 



where the Wilson-line operators (0, -F ei i,0)" and 
(0, G e i,0)" are the operators © made from A and B 
fields, respectively. As discussed in 0, 0, |23|, 0] , the 
slope of Wilson lines is determined by the "rapidity di- 
vide" vector e Vl = pi + e~ Vl p2- (From the wiewpoint of 
A fields, the slope e\ can be replaced by p% with power 
accuracy so we recover the generating functional © with 
(0,G H ,0)=U u -U 2 i). 





FIG. 4: Effective action as a scattering of two shock waves. 

III. SCATTERING OF OCD SHOCK WAVES 

A. Efffective action as a shock-wave scattering 
amplitude 

In this section we define the scattering of the shock 
waves using the rapidity factorization developed above. 
Applying the factorization formula (|13f) two times, one 
gets (see Fig. 0J: 



DA J( P a)J(pb)J(-p' a )J(-p' b ) e lS ^ 



(14) 



iS(B) 



J DAJ{p A )J{-p' A )e iS ^ J DB J( PB )J(-p' B ) e l 
DCexp \iS{C) + iJ d 2 z ± {{0 : A eil ,0} a z {0,C eit ,0} a z 

+ (0,^,0)2(0,^,0);} 

where the slope is e\ = P i + e~ m p2 for the [...] Wilson 
lines and e2 — pi + e~ m p2 for the (...) ones. 

The functional integral over the central range of rapid- 
ity rji > 7] > r\2 is determined by the integral over G field 
with the sources 

(0,A ei i,0) z = [0,coe 1 ] z (id i + gAi{ooe x + z±))[ooe 1 , 0] z 

- [0, -ooeA z (idi + gAi(-ooei + zx))[—ooe 1 ,0] z 
(0,B e2l ,0) z = [0,coe 2 ] z (idi + gBi(ooe 2 + z_i_))[ooe 2 , 0} z 

- [0, -ooe 2 ] z (idi + gBi(-ooe 2 + zj_))[-ooe 2 , 0] z (15) 

made from "external" A and B fields. Since -Aj(±oo) 
is a pure gauge these sources can be represented as a 
difference of a pure-gauge fields (0, A e i,0) z = (V™ — 
Vf) z and (0,B e2l ,0) 2 = (U? -Uf) z where 

Vi,2(z±) = [0, ±ooei] z [±ooei + z±, ±ooei + ooe^] 
Uifl{z±) = [0, ±ooe2] z [±ooe2 + z±, ±ooe2 + ooejj (16) 

Since there is no field strength F^ v at infinite time the 
direction of e± does not matter. 

The result of the integration over the G field is an 
effective action for the r\\> r\> r\2 interval of rapidity 

e iS ff (Vi,V 2 ,£/i,!72;7)i— 172) ^7) 

= ^Gexp [iS(C) + tJd 2 z ± {(V? 1 - V 2 M ) 2 [0,G eiI ,0]: 

+ (Uf -u?uo,c e2i ,o) a z }~ . 
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One can interpret Eq. I|17|) as an effective ac- 
tion for scattering of two QCD shock waves denned 
by the sources (|lf>|> . Note that the effective action 
iS ef f(y 1 ,V 2 ,Ui 1 U 2 ;r]i ~ m) defined by Eq. JT7|) is in- 
variant under the redundant gauge transformations 10 



U\(2){xi_) -> u 1(2 ){x±)n(x±), 



(18) 



since this transformation can be absorbed by a gauge 
rotation of the C fields 



C M - ^ (x ± , In — )C^(x x , In — ) 

+ -n\x^\n—)dMx^M—) 

g X, x. 



(19) 



where Q(x_i_,ln^) is an arbitrary SU3 matrix satisfy- 
ing the conditions fl' (x±,r)i) — 0(xj_) and £1^ (ecj^ , 172) = 
<)'(.<• !• 

With a power accuracy 0(m 2 /s), we can replace e\ by 
pi and e 2 by p 2 : 



(Vr-Vn^O.Ci.O]" 



(20) 

DCexp {iS(C) + i Jd 2 z ± [(Vf - V™) z [0, C. h 0]° 

+ (wf-z^co.Oi.o);]} 

The saddle point of the functional integral l|2l)|) is de- 
termined by the classical equations 

+ (Uf -U?) z (0,C*i,0):]} = (21) 

At present it is not known how to solve this equations 
(for the numerical approach see 0] ) . In the next section 
we will develop a "perturbation theory" in powers of the 
parameter [U, V] ~ g 2 [Ui,Vj]. Note that the conven- 
tional perturbation theory corresponds to the case when 
Ui, Vi ~ 1 while the semiclassical QCD is relevant when 
the fields are large (Ui and/or V,; ~ i). 



_6_ 

6C, 



B. Expansion in commutators of Wilson lines 

The effective action is defined by the functional integral 
(|20|l (hereafter we switch back to the usual notation 
for the integration variable and for the field strength) 

giSeff (Vi,Vi,l7i,l72;i7i— 772) 

= jDAcxp (iS(A) + 1 fd 2 z x {(Vf - VrUO, F.i, 0} a z 

+ (U?-U?U0,F*i,0)z}) (22) 

Taken separately, the sources ~ Ui create a shock wave 
Uu6(x*) + U 2 i0{— x*) and those ~ Vi create Vu6(x,) + 



V 2 i6{—x m ) In QED, the two sources Ui and Vi do not 
interact (in the leading order in a) so the sum of the two 
shock waves 



f(0) 



: U u 0(x.) + U 2l 8(-x*) + V u e(x.) + V 2l 8(-x.) 



a: 



(0) _ 1(0) _ 



a';" = 



(23) 



is a classical solution to the set of equations lf?T|l . In 
QCD, the interaction between these two sources is de- 
scribed by the commutator g[Ui,Vk] (the coupling con- 
stant g corresponds to the three-gluon vertex). The 
straightforward approach is to take the trial configura- 
tion in the form of a sum of the two shock waves and ex- 
pand the "deviation" of the full QCD solution from the 
QED-type ansatz l|23|) in powers of commutators [U, V] . 
This is done rigorously in |l6j | and the relevant formulas 
arc presented in the Appendix. Here wc will use a slightly 
different zero-order approximation (cf. []]) which leads 
to same results in a more streamlined way at a price 
of some uncertainties (like 6(0)) which, however, do not 
contribute to the effective action in the leading order. 

Let us consider the behavior of the solution of the YM 
equations at, say, xq — > 00, X3 fixed (in the forward quad- 
rant of the space). Since there is no field strength at 
t — > 00, the field must be a pure gauge. As demon- 
strated in Ref. 0], this pure-gauge field has the form 
of a sum of the shock waves plus a correction propor- 
tional to their commutator. Technically, for a pair of 
pure gauge fields Ui(x±) and Vi(x±) we define Wi(x±) = 
Ui(x±) + Vi(x±) + gEi(x±;U,V) as a pure gauge field 
satisfying the equation (idi + g[Ui + Vi, )E % = 0. In the 
first order in [U, V] this field has the form 



E?(U,V) = -(x ± \U^ T U^ + vK-V* 
PI P± 

P k 1 oft 



\ ab [U h V k ] b -z^fc)+0([C/,y] 2 ) (24) 

p± 

where [Ui,V k ] a = 2Tit a [U t ,V k ]. The second, [U, V] 2 , 
term of the expansion (|24|) can be found in |lq | but we 
do not need it with our accuracy. 

Throughout the paper, we use Schwinger notations for 
the propagator in the external field (x\^\y). For the 
bare propagator it reduces to (x\^\y) and for the two- 
dimensional propagator in the transverse space we use 
the notation (x±\ ^-\y±) where p\ — ~PiP % ■ Also, |/) 

denotes Jd 2 z±f(z±)\z±) and later we will use the nota- 
tion \0,f)=Jd 2 z x f(z ± )\0,z ± ). 

The zero-order approximation for the solution of the 
classical equations for the functional integral (12211 can 
be taken as a superposition of pure gauge fields in the 
forward, backward, left, and right quadrants of the space 
(see Fig. |5J|: 



A[ 0) 



a! 0) = 

Wj.(x ± )0(x*)0(x.) 



(25) 



Wi(x±)6(-x.)6(x.) 



W r (xa_)6(x x )6(-x.) + W l B (x ± )8(-x*)8(-x.) 
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FIG. 5: Pure-gauge ansatz as a zero-order approximation to 
the classical field created by the collision of shock waves. 



where 

W i F =Ui+V{ + E i F , W l L =U l 2 + V\+E F 
W i R = U{ + V i 2 +E F , W B =Ui+Vi+E F (26) 

and E F {Ux , Vx) , (U 2 , Vx) , E R {U X , V 2 ) , and E B (U 2 , V 2 ) 
are given by Eq. Ij24(l . For the trial configuration (|25|) 

PS = 6(^){9(x.)(W F - W L ) + 9(-x.)(W R - W B )} 

2r 

K = 6(^){6(x*)(W F - W R ) + e(-x*)(Wh - W B )} 



(27) 



so 



D l F, 



\i = 5{-xM[e{x.){d l -i[A\){W Fl - W Ll ) 
s V 



D l F~i = 



+ e{-x.){d l -i[A\)(w Rl -w Bl ) 

8(-x.)(e(x*W - i[A\)(W Fi - W Ri ) 
s V 

+ B(-x.)(8 i -i[A i ,)(W Li -W B i) (28) 



and 



D*F.i = D.F mi (29) 

= 5(-x*)5(-x.){W Fl -W R ,-W Lt + W Bl ) 
s s 

2 2 

= S(~x*)6(~x.)(E Fi - E Rl - E Li + E Bi ) 
s s 

Next, one shifts A — > A+A^ in the functional integral 
(J22J) and obt ams 



(30) 



/" cxp + i J d i z{^D lxv A v + rM„)}. 



Here 



.{(Vx - V 2 )?(W> - w£ + - W^) 4 



+ {Ux - u 2 )t(w F + w L - w R - w B y a - -{w F 
-w L + w R - w B y a (w F + w L -w R - w b ) c i (3i) 



is a sum of the action and source contributions due to 
the trial configuration — D 2 (A)g lJLl/ — 2%F ViV is 

the inverse propagator in the background- Feynman gauge 
|32| and T M is the linear term for our trial configuration: 



Ti = 25(-x»)8{x.){W F -Wl-W R + W B ) 
= 26(~x*)6(x.)(E F -Eh-E R + E B ) 



I 



s 

l6(^x,)(6(x*)[Vxi - V 2i ,E F + E R ] 
+ 9(-x*)[Vxi-V 2i ,Ei + E B ] 
T. = -U{-x*)(6{x.)[Uxi-U 2i ,E F + E L ] 
+ 6{-x*)pxi-U 2i ,E R + E B t 



(32) 



The first line in this equation follows directly from Eq. 
(|29|l while the two last lines are obtained by adding Eqs. 
(|28|l and the corresponding first derivatives of the sources 

ED 



6 A, 



d 2 z±(Uf -W 2 a %(0,F„,0)^ 



A,=0 



JaJ^ 



= £(£*){ - W°OP2 + Xl_),)(Uxi - U 2i ) 
+ 9(-x.)(d i - i[A\-oop 2 + xx),)(Uxi - U 2i )] 

"(Vr-V 2 ai ) 2 [O.F.i.O]" 

A.=0 

= 5(x.){e{x*)(d l - i[A l (^ Pl + x± ), )(Vxi - V 2l ) 
+ e{-x*W - »[^(-oopi + x ± ),)(Vxi - V 2l )} (33) 
We get 

^-{Sqcd + J d 2 z ± (U?-Ur) z (0,F„,0r z } 

= - 8(x,){d(x.)(8 i - i[A l (oop 2 + x x ),){E Fi - E Li ) 
+ 9{-x.W - i[A\-^ P2 + x ± ),)(E Ri - E Bi )}, 

^-{^qcd + J d 2 z ± [(Vf - V§%[0, F. h 0]«} ^ 

: -8(x m ){o(x*)(8 l -i[A i (oop 1 +x±),)(E Fi -E IU ) 

+ 9(-x*)(d l - i[A l {-oopx + x ± ),)(E Li - E Bl )} (34) 

Using 0(0) = | so that A i {oop 2 +x x ) = \ {W m + W Bl ), 
A l (-oop 2 + x ± ) = \{W Li + W Bl ), A^oopx + x ± ) = 
\(W Fl + W m ), A^-oopx + x±) = \(W Ll + W m ), and 
the condition (idi + [Wi, )E l — one easily obtains Eq. 

Expansion in powers of T in the functional integral 
(|30|l yields the set of Feynman diagrams in the external 
fields l|23|l with the sources (|32|l . The parameter of the 
expansion is g 2 [U t , Vj] (~ [U, V], see Eq. (J5J). 
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IV. THE EFFECTIVE ACTION 
A. The effective action in the lowest order 

The effective action (|22[1 in the first nontrivial order 
in [U, V] is given by the integration of linear terms i|32|) 
with the Green functions in the external field (|25|l 

iS eS (U,V) = 

-\ [ d i zd i z l T a (^ a (z)A^(z'))T b u (z l ) (35) 



It is easy to see that the term ~ T»T. is ~ [U, V} 3 so the 
leading contribution ~ [U, V] 2 comes from the product 
of two Tj's which has the form 



2 

where 



d 2 z ± d 2 z' ± mz ± )(0,z ± \ 



1 



U 



2(E i F - 
2{W l F 



E' T 



p2 



El 



-\0,z> ± y»L b \z> ± ) 



E l B ) 



(36) 



(37) 



is actually the transverse part of the Lipatov vertex of 
the gluon emission by the scattering of two shock waves 
in the first order in [U, V] (see Appendix). As we shall see 
below, the main logarithmic contribution to the integral 
(|36[1 comes from the region z± — > z'j_ where one can re- 
place the propagator in the background field by the bare 
propagator. One obtains 



dad(3 
~8^~ 



(0,i?| 



a /3 s 



p\ 



10, V 



(38) 



The integral l|38|l is formally divegrent. Within the 
LLA accuracy, one can cut the integration off at the 



width of the shock waves A 
and obtain: 



»W2 



dad/3 
a s Ar) 



-i(a\+(3p) 



(0,i?| 



a/3s - 



Pi 



it 



|0,L m ) 



d 2 z ± Lf{z 1 _)L la {z 1 _) 



(39) 



where Arj = rji — 772 is our rapidity interval. 

In addition, within the LLA approximation the zero- 
order term (|31J) can be simplified to 



(40) 



S = d 2 z ± {Vi-V 2 ) m {U 1 -U 2 ) a i 



Indeed, it is easy to see that in the r.h.s. of Eq. I|31|l the 
terms ~ [U, V] cancel while the the terms ~ [U, V] 2 are 
not multiplied by Arj so they can be omitted in the LLA 
(the Eq. is ~ [U,V] 2 Ai]). It is worth noting that 

Eq. I|4(J|) is the usual light-cone lattice action in the 
limit when transverse size of the plaquette vanishes and 



the longitudinal increases to infinity. Thus, the effective 
action in the first order can be represented as 

S eS (U,V) = jd 2 z^[{V l -V 2 ) a l {U 1 -U 2 ) a l 

- J d 2 z ± L?(z ± )L™(z ± )} (41) 

We shall see below that L; is the Lipatov vertex of the 
gluon emission by the scattering of two shock waves in 
the first order in [U, V] . Note that S c s given by Eq. |gT} 
is invariant with respect to rotation of the sources 



(42) 



For the first term in the r.h.s. of Eq. I|41|) it is trivial 
while for the second it follows from the gauge-invariant 
form discussed in the next section, see Eq. I|56|) . 

For future applications we will rewrite the effective ac- 
tion (|41|) as a Gaussian integration over the auxiliary field 
A coupled to Lipatov vertex i|A9(l : 

e iS eS (u,V) _ e iJ d 2 z ± (Vi-v 2 y a (Ui-u 2 )? (43) 
x J DXexp { - a s A V J d 2 z±(XfX ai - L"A ai ) | 



B. Nonlinear evolution equation from the effective 
action 

Let us prove now that the effective action (|4^|) agrees 
with the non-linear evolution equation. To find the evolu- 
tion of the dipole U X U^, we need to consider the effective 
action for the weak source V. From eq. (|24|l one sees 
that at small gVi ~ diV 

Lf(x ± ) = -2{x\ut P -^U t - U^U^iVr - V 2 f k ) 



(44) 



and Eq. I|22[l can be rewritten as 

fDAexp {iS(A) + 1 J d 2 z ± [(Vr - V 2 m ) z [0, F.i,0] a z 
+ (U?-U?) z (0,F*i,0):]} 



2 Jz\ 

Here 

Ui = e d i Ui, 



2 as A V -^-(u 2 yul) 
U 2 = e d i U 2 



so that 



U a u = U? l + 2a s A v (Ut^ r U 1 r»\ b k 



did k , 
. did k _ 



U« i + 2a s A v (U^U 2 r h \ k 



(45) 

(46) 
(47) 

(48) 
(49) 
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Wilson loop is defined as follows 




FIG. 6: Gauge-invariant form of the effective action. 



where we need only the first term in expansion in A; in 
A. 

To find the evolution of the color dipole J7J) we should 
expand Eq. in powers of Vu — V 2 i and use the for- 
mula 



[f^oopiyau +oopi,y_L + oopi][ocpi,0]^ 2 
= p e i 9n dz i[°> F '*> (S \* 

which results in 



(50) 



DX e sA vI ^ z± xtx" u f ly U lx ulU2y (51) 



Performing the Gaussian integration over A one obtains 
after some algebra 



tr{U lx UlU 2y Ul y } 



(52) 



(z - y)± 



= tr{U lx uLu 2y Ul v } + ^ /«<V {x _ z)l{z _ y)2± 

x (tr{U lx Ul U 2z Ul }tr {U lz Ul U 2y U\ y } 

-N c tT{Ui x UlU 2y Ul y }) 



which is the non-linear evolution equation ^26, 27] for the 
Wilson-line operator U x = Ui x U 2x = [ooe, — ooe\ x . 



C. Gauge-invariant representation of the 
first-order effective action L;L l 



Our expression for the A77 term in the effective action, 
proportional to the square of the Lipatov vertex (Wp — 
W L — W l R + W l B ) a , was obtained in the axial-type gauges. 
It can be rewritten it in the gauge-invariant "diamond" 
form of trace of four Wilson lines at av* = ±00 (see Fig. 

as suggested in a recent paper 0]. The "diamond" 



0(£_l) = tr{[-oopi,F. i ,oop 1 ] oop2 [oop 2 ,F* i ,-oop 2 ] oopi 

X [00pi,-00pi]_ OO p 2 [-00p 2 ,00p 2 ]-oop 1 

+ tr[-oopi, oopi] oop2 [oop 2 , F H , -oop 2 ]oo Pl 
x [oopi, F,i, -oopi]_oo P2 [-oop 2 , OOp 2 ]-oopi 
+ tr[— oopi, oopi]oop 2 [oop 2 , — oop 2 ] 
x [oopi, F,i, oopi]- oopa [-oopa, F*i, ooP2]-oo Pl 

[oop 2 , -oop 2 ] oopi 
x [oopi,-ocpi}- oop2 [-ocp 2 ,F*i,ocp 2 ]-oop 1 } (53) 

where the transverse arguments in all Wilson lines are 
xx- Next, define this "diamond" as a function of the 
sources 

0(U lt U 2 ,V 1 ,V 2 ) = (54) 
M-'JdA 0(A) exp (iS(A) + i J rf 2 z ± {(vr - Vf) z 

x [0, F eiU 0]° + (Uf ~ Uf) z (0, F e2U 0)°}) 

(In the Appendix we demonstrate that the trace of four 
Wilson lines 

tr{[— ooei, ooei]ooe a [°° e 2, -ooe 2 ]oo ei 

x [ooei, -ooeij-coezC- 00 ^, ooe 2 ]-oo ei } = 1 (55) 

is trivial in the leading order.) 

Note that ()(U,V) is invariant with respect to the ro- 
tation of all sources by one gauge matrix il(x±) 

0(u 1 n,u 2 n,v 1 n,v 2 n) = o(u 1 ,u 2 ,v 1 ,v 2 ) (56) 

since it can be absorbed by gauge transformation of the 
fields A^ — > Q'A^Q, + ^f^9 M Q in the functional integral 
Eq. 

Now we can prove that the square of Lipatov vertex 
can be expressed as the "diamond" Wilson loop: 



1 



L1(U,V)L ai {U,V) = 0(U,V) 



(57) 



Indeed, it is easy to see that for the trial configuration 
the Eq. reduces to 

(W Fi - W Li ) a (W l F - W l R ) a + (W Ri - W Fl ) a 

x (W R - W B ) a + (W Bt - Wfli)°(Wfl - wD a 

+ (W Lj - W Bi ) a {W L - W l F ) a = 

[W Fl - W Lt -w r + w Bl ) a (w l F -W l L -W R + W B ) a 

which coincide with the l.h.s. of the Eq. I|57|) . In the 
covariant-type gauges where Ai — + as xu — > oo the r.h.s 
of the Eq. (|53ll can be rewritten as 

d{ Mi d i M 2 Afj m| + MxdiM^MlMl 

Mi M 2 di m\ di m\ + d l MiM 2 Mld l M\ (58) 

where Mi = [-oopx, oop 1 ] oop2 , M 2 = [oop 2 , -oop 2 \ oopi , 
M 3 = [oop^-oopij-oo^, and M 4 = [-oop 2 , oop 2 ]^ oopi . 
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Eq. (|58|l is the expression obtained recently in |l7l in 
the framework of the Hamiltonian approach (see also |l9l 
[2£j ) ■ The corresponding form of our effective action is the 
following: 

hf(U,V)L ai (U,V) 

= di{W L W F )d\W F W R )W R W ] B W B wl 

+ WLWld^WpW^diiWRW^WsWl 

+ WLWlWpW^iWRWl^iWBWl) 

+ ^(WLW^WpWlWRWld^WgWl) (59) 

The Eq. (|57|l links the representation in terms of the 
effective degrees of freedom (Wilson lines in our case) 
with the representation in terms of gluons of the under- 
lying Yang-Mills theory via Eq. (|54|l . The remarkable 
feature of the gauge-invariant form 1)57(1 is its universal- 
ity - if one writes the effective action in terms of some 
other degrees of freedom (say, reggeized gluons (2^) one 
should recover Eq. I|57(l once these new effective degrees 
of freedom are expressed in terms of gluons. 




FIG. 7: Effective action factorized in n functional integrals. 



V. FUNCTIONAL INTEGRAL OVER THE 
DYNAMICAL WILSON LINES 



A. Effective action as the integral over the Wilson 

lines 



In this section we will rewrite the functional integral for 
the effective action (122ft in terms of Wilson-line variables. 
To this end, let us use the factorization formula ((13(1 n 
times as shown in Fig. [7| The effective action factorizes 
then into a product of n independent functional integrals 
over the gluon holds labeled by index k: 



JS e![ (U,V;v) - 



(60) 



+ s(A n+1 ) + (vr" 1 ' 1 - vj +1,l )(wr - W 2 m ) 

+ S(A n ) + ... + (Vf - Vf){U l u - Z4) + S(Ar) 
+ {V\ l -Vl l ){Uu-U u )} , 



where the integrals over xj_ and summation over the color 
indices are implied. As usually, U k ' z = ^U^ k d % U k and 
y k,i = iyikgiyk where 

VlCqOs-l) = Pe l9 ^ x ' du <-^ Ak ^ un '" 1+Xl - ) (61) 

and the vectors n k are ordered in rapidity: r) Q > rj n > 
r\ n -\ ■ ■ ■ r]2 > Vi > Vo- To disentangle integrations over 



different A , we rewrite exp^l/^ 



fc+M 



V. 



U^ 1 )] at each "rapidity divide" ijk as an integral over 



the auxiliary group variables V^ 1 and U k 2 using the 
formula 



(62) 



x / DV(x ± )DU(x ± )e^ ' dx ± v,u' +l fdx^w-i f dx ± v,w 



(where Vi = V^diV and U l = U^diU). The determi- 
nant gives the perturbative non-logarithmic corrections 
of the same order as the corrections to the factorization 
formula 1(13(1 . In the LLA they can be ignored, conse- 
quently, we obtain 



iSet{{Ui,U2, Vi,Vb,J7i —772) 



j DA 1 DA 2 . . . DA n+1 exp i { (Vj - V l 2 ) (W™ +1 - U^ 1 ) = j 



(63) 



Ul+lDA k Iil =0 DU k DU k DV^DV 2 k 



x exp 1 



{(Vi - V 2 ){Ul +1 - U^ 1 ) + S(A n+1 ) + (V£ +1 

+ (vi - vDiuf ~ uf) - (vf - vfwi - ul) 

+ (V| - V 2 u ){ur - Uf ) + S(A 2 ) + (V 2 U - Vl){U? 

- uf) - (% - - u?) + (vh - %) 

x {U\ l - Uf) + SiA 1 ) + (Vl l - Vf)(U lt - U 2i )}. 

Now we can integrate over the gluon fields A k ■ Using the 
results of the previous Section, we get 



DA k exp{z(V 1 fc I - V k 2l ){U k S - U k f) + zS(A k ) 



+ *(v k l ~v k l )(u k - 1 ' l -u k - 1 ^)} 



(64) 



9 



where at sufficiently small Arj 



i?^L i {V\U k - 1 )L\V h ,U k -i) (65) 



Performing the integrations over A k we get 



x exp 1 



{(vj-vM" 



z4 



- (v 2i - v 2i )(z4 - z4) + (v?, - v 2 )(z^ - 



^L 2 (f 2 , ^)Ar? - (Vi - v£)(ft" - W a ") 

+ (Vi, - V 2 \)(^ - f/ 2 i) - ^L 2 (^\ C/)A^}. (66) 

In the continuum limit n — > 00 we obtain the following 
functional integral for the effective action 

e iS a tf(U 1 (x),U 2 (x),V 1 (x)y 2 (x); m ~ V2 ) 



(67) 



n 



2W J (i,l)W ) (x,')) 



x exp 



d 2 x (itV^Cs) - V^)WfM -U?(x,ri)} 



x — [Uf{x,<n)-U¥{ X ,ri)} 



^-Lt(V(x,r,),U(^v))L ai (V(x,rj),UM)}) 



where we displayed the color indices explicitly and re- 
moved the hat from the notation of the integration vari- 
ables. This looks like the functional integral over the 
canonical coordinates U and canonical momenta V with 
the (non-local) Hamiltonian L 2 (V,U). The rapidity 77 
serves as the time variable for this system. The above 
representation of the effective action as an integral over 
the dynamical Wilson-line variables is the main result of 
this paper. 

Note that the LiU term in the exponent in (|67|l is 
invariant under the rotations (|18|l 

Uj(x,ri) -> [7j(x,?7)n(a;,?7), Vj(x,ri) -> Vj (x, rj)£l(x, rj) 

(68) 

(see. Eq. <|56|0 . but the term ~ U-^-V preserves only the 
77-independent symmetry 

Uj(x,r)) Uj(x,r))Q,(x), Vj(x,T)) -> Vj(x,r])Cl(x). 

(69) 

This probably means that the term ^ U-^V should be 
adjusted by a ~ [U, V] 2 correction (not important in the 
LLA) so that the full symmetry (|68|l is restored. 




FIG. 8: Feynman diagram for a typical [U, V] 3 term. 



The idea how to use the factorization formula to 
rewrite the functional integral in terms of Wilson lines 
was formulated in Ref. Q] where the first-order effective 
action was obtained (the expression in terms of square of 
Lipatov vertex is given in [l6j). However, the additional 
redundant gauge symmetry (|18|) was fixed by the require- 
ment that there is no field at t — > —00 which correspond 
to the choice U2 = and V 2 = for the two colliding 
shock waves. In this case, one obtains the functional in- 
tegral in terms of only two variables, U and V, at aprice 
of a more complicated form of the effective action [lj . 

It should be noted that L 2 (U, V) is only the first term 
of the expansion of the true high-energy effective action 
K(U,V) in powers of [U, V]. An example of the next- 
order, ~ [U, V} 3 , contribution to K(V, U) which is miss- 
ing in the effective action ijoTjl is presented in Ref Q, see 
Fig. 



B. Functional integral for the non-linear evolution 

It is instructive do demonstrate that the functional in- 
tegral l|67|) reproduces the non-linear evolution in the case 
of one small source. Basically, we recast the arguments 
of the Sect. IIV Bl in the language of functional integrals. 

First, note that at small V the functional integral over 
V is Gaussian (see the Eq. (|44l) i. It is convenient to in- 
troduce the "gaussian noise" associated with the Lipatov 
vertex and rewrite the functional integral (|67|l as: 



JS B f { (U 1 (x),U 2 (x),V 1 (x),V 2 (x);r !l -r l 2) 



(70) 



D\?(X,T)) 



= fn j=h2 DV ] (x,7 1 )DU J (x,r 1 ) 

J C/j(x, 77 2)=!7 i (x) 

x exp Jd 2 x{i[V^(x)~V^(x)}[Ur(x,v)~H2 l (^v)} 
rvo ( 



- i[VUx,v) - V^{x,rf)\ — [Uf{x,ri) - Uf(x,r,)] 

- 2a s X ai (x,v)(x\UI^U 1 - (1 <- 2)r b (V? - 

(When convenient, we use the notation (...)'' = (...)(rj) 
for brevity). The integral over V gives the (5-function of 
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the form 



6(di(— [Uf(x,r,)-Uf(x,ri)} 

k k 

- 2ia s (x\U[ — —Ui - (/^— 2-(/ 2 | X k ') 

p ± 



Pa 



which restricts U in a following way 
d 



th] 



[U?(x,ri)-U?(x,ri)] 



= 2i as {x\V\ VJ ^Vx-V\ V -^V2\ ah \t) (71) 
Px Px 

It is convenient to rewrite Eq. (|T1|) in the integral form 
(cf. Eq. H7J): 

U l {x,r)) = Te p ± Ui(xx,V2) 

(72) 

where T means ordering in rapidity (= our "time"). The 
remaining integral over A is gaussian with the "propaga- 
tor" 

W{x ± ,rj)\Uy ± , */)> = 9ijS ab ^-6(xx yx)8(r 1 rf). 



2a„ 



The evolution of the dipole can be represented as 



(73) 



The factor 1/2 in the Eq. O comes from 9(0) = 1/2. 
To avoid this uncertainty, one should first calculate the 
correlations in A and then differentiate with respect to 
rapidity (cf. Ref. ^ ) 

§jfdr)'dr," (A™ (x x , rf)X b j (yx , rf ) 



2a 



-S(x - y)x9ijf(v)g(v) 



|- fdrl f d^iAt&x^^V^fiv'W) 



1 

4a. 



■6(x - y)xgijf(v)f(v) 



The first line in the above equation should be used to 
make contractions between different T and T in Eq. I|74|) 
while the second line takes care of the contractions within 
same T or T. It is easy to check that the result is con- 
sistent with taking 0(0) = 1/2 in the Eq. 1(751) . 

Similarly one can demonstrate that all the hierarchy 
of the evolution equations for Wilson lines j2|| |3(j (= 
JIMWLK equation [^H]) is reproduced. 



C. Classical equations for the Wilson-line 
functional integral 



12 ' ' • 

i„ - -2iast a f n d V '(x\j£-U'/\ ab \''/) 

2ia 3 t a n dr 1 '(y\li-U rl ' \ ah \ b l n ' ) „ 

x Te Jv2 ^ k Uf y (74) 

where T denotes the inverse rapidity ordering. Taking 
the derivative with respect to r\ we get 



d 



— TP 'Tf ' TP HI 

d<n v 2y 



DX 



-a a Qd v fd 2 z ± X^\^ 



x u 2x u 2y 



x 2ia s (up y [t a (x\^\ ab Xf) - x - y]Ul 
v Px 

- U%Ulu£[nx\^U^Xt>) - x ~ y]Ul) 

Using the contraction 

(\%{z,n)U{x,ri)) = -\i(z\U\r,)^\x) ab t b U(x,r,) (75) 

(Xl(z, V )U\x,r,)) = -h( z \uH v )^\xy b U\x, V )t b , 
z p ^ 

one gets the non-linear evolution equation l|52|l after some 
algebra. 



As we discussed above, the characteristic fields in the 
functional integral are large but the coupling constant 
q s (Q s ) is small due to the saturation. In this case, we 
can try to calculate the functional integral l|t)7|) semiclas- 
sically. Using the approximate formula 

8W(U, V) ~ -(W^^-W) ab (SU b i + 5V b ) (76) 
Px 

we get the classical equations for the functional integral 
(|67[) in the form 



(d i -igVi) ab (U u -l(2i) b 



(77) 



= 2ia s (d i - igVi) ab (Wl^-W F - W t P -^W L ) bc E c r 

Px Px 

(d i -igVi) ab (U u -l{ 2i ) b 



= 2ta s (d l - ig V l 2 ) ab (W R P -fw R - W B P -^W B ) bc E c r 

Px 



(d l - tgU{) ab (V u - V 2l ) b 



= -2ia s (d l - igU{) ab (Wl V ^-W F - wIP^Wr^E], 

Px Px 

(d l ~igW 2 ) ab (V u -V 2l ) b 

= -2%a„{& - igVl) ab (wi P -^W L - wl v y~W B ) bc E] 

Px Px 

with the initial conditions 

U(rj) = U at r] = r? 2 , V(rj) = V at r) = rji. (78) 



Px 



p l p> . 
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At small Vi these equations reduce to (cf. Eq. 171(1 ') 

{Uu-U^Y = 2za s ([/ 1 t ^C/ 1 - ul^-U 2 ) ab E) 
PI PI J 

Vu-V 2i = 0([U 1 V] 2 ) (79) 
while in the opposite case of small Hi they are 

(Vu - V2i) a = -2ia s (V^V! - V 2 f ^V 2 ) ab E b 

P± PI 

Uu-U 2l = 0([U,V] 2 ) (80) 

It is instructive to rewrite the equations (|79|l and (|80|l in 

terms of W's. 

Wk - W a Bi = 2ia s {W^-W R - wl^W L ) ab E c 3 
Pi Pi 

Wm - W a Ll = 2ia s (Wl^W F - wl^W B ) ab E^ 
P± PI 

(81) 

From the viewpoint of the functional integral (|67fl the 
W's are the (non-local) functions of U and V variables 
given in the first order by Eq. I|24[l . It would be very 
interesting to rewrite the Eq. IL'l'li of the W variables 
themselves, that is, to construct the functional integral 
over the W variables with a saddle-point equations given 
by the Eq. (gTJ. 




FIG. 9: Wilson-line structure of the effective action . 




FIG. 10: Typical Feynman diagrams included in the effective 
action. 

pomeron loops, see Fig. ^3 I n the dipole language, it de- 
scribes both multiplication and recombination of dipoles 
(see the discussion in [l8|, W$ ) ■ In conclusion I would like 
to emphasize that the effective action (|67(1 summarizes 
all present knowledge about the high-energy evolution of 
Wilson lines in a way symmetric with respect to projec- 
tile and target and hence it may serve as a starting point 
for future analysis of high-energy scattering in QCD. 



VI. CONCLUSION 

As mentioned in the Introduction, the popular idea of 
how to solve QCD at high energies is to reformulate it 
in terms of the relevant high-energy degrees of freedom - 
Wilson lines. The functional integral (|67|l gives an exam- 
ple of such 2 + 1 theory where 2 stands for the transverse 
coordinates and 1 for d the rapidity serving as a time vari- 
able. The structure of the effective action is presented in 
Fig. |5J Note that the two terms in the exponent in the 
effective action, shown in Fig. |SJ are both local in x± 
but differ with respect to the longitudinal coordinates: 
the first (kinetic) term is made from the Wilson lines lo- 
cated at x + = or .x_ = while the second term is made 
from the Wilson lines at x± — ±oo. Unfortunately, the 
transition between these Wilson lines is nonlocal in x± 
(see Eq. (I24|) ') and so the resulting effective action is a 
non-local function of the dynamical variables U and V. 

In should be emphasized that Eq. i|tj7fl is only a model 
- the genuine effective action for the 2 + 1 high-energy 
theory of Wilson lines must include all the contribu- 
tions ~ [U, V] n (as we mentioned above, an example of 
a [U, V] 3 term which is missing in Eq. I)22|) is presented 
in Fig. |SJ). However, this model is correct in the case of 
weak projectile fields and strong target fields, and vice 
versa. In terms of Feynman diagrams, the effective action 
(|67|l includes both "up" and "down" fan ladders and the 
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APPENDIX A: CLASSICAL FIELDS IN THE 
FIRST ORDER IN [U, V] 

Following Ref. |l6j . we take the zero-order approxi- 
mation in the form of the sum of the two shock waves 

Af ] = UiB{-x m ) + Vi6{-x.), A { ° ] = il 0) = 0(A1) 

We will expand the "deviation" of the full QCD solution 
from the QED-type ansatz I|A1|I in powers of commuta- 
tors [U, V]. To carry this out, we shift A — > A + A^ in 
the functional integral l|22|) and obtain 

j DA cxp{i j d 4 z(±A»D^A» + 5 T"A M )}.(A2) 
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Here = D 2 {A)g ilv — 2iF^ v is the inverse propagator 
in the background-Feynman gauge and is the linear 
term for the trial configuration I|A1(I . Since the only non- 
zero component of the field strength for the ansatz (|A1|) 
is 



F, 



(0) 



ik 



-i[U u ,V lk ]d(F) - i[U u ,V 2k ]9(R) 



(A3) 



i[U 2l , V lk ]9(L) - i[U 2i , V 2k }9(B) - (t «-> k) 



the linear term T M = DPF$ is 



T* = T. = 0, 

r = ( i a fc + .g[c/ lfc + T/ lfe ,)([f/i,y 1 fc ] 



(A4) 
z <-> fc)0(F) 



i fc)0(L) 
i «-» fc)0(B) 



+ (id fc +.g[l/ 2fe + U 2fe ,)([(7 2 \Vf 

+ ^ fc + . g [c/ 2fe + y 2fe ,)([^,y 2 fc 



where 61(f) = 0(z*)0(z.), 0(i?) = 0(z*)0(-z.), 0(L) = 
9(-z*)6(z.), and 0(B) = 9{-z*)9{-z.) 

Expanding in powers of T in the functional integral 
(|A2I) one gets the set of Feynman diagrams in the external 
fields H23[) with the sources l|A4l) . The parameter of the 
expansion is g 2 [Ui, Vj] (~ [U, V], see Eq. (J5J). 

The general formula for the classical solution in the 
first order in [U, V] has the form 

4 1)a W = igj d A z(Al{x)A b »{z)) A T b u {z) (A5) 

The Green functions in the background of the Eq. 123|) 
field can be approximated by cluster expansion 

(A^x)A»(z)) A 

= {A M {x)A v {z))u + (A^x)A»(z)) v 

- (A^x)A v (z)) + O(\U,V\) (A6) 

where {A^(x)A v (z))o is the perturbative propagator and 
(Al{x)Al{y)) v x * >0 ^ <0 / dz8(-z*) {ul(x\^—\z) 

p J s I + l€ 

2ag llv U\ U\ + — (p 2v d tl {U l Ul) z + /i <-> i/) 



/ jr + ze / 



(A7) 



is the propagator in the background of the shock wave U 
(the propagator in the V background is obtained by the 
replacement U <-> V, pi Pi ■ 

Substituting Eq. (|A7|I and l|A4J) into the above equa- 
tion, one obtains (the details of the calculations can be 
found in Ref. |16| and here we present only the the final 
set of gauge fields): 

A" - 9(F){w^(x^) - gt a (W F -^—W F ) ab L F b (x)} 
+ (F <-> L) + (F <-> R) + (F <-> B) (A8) 



where L F = L l L — L R = U B = 2E l and 



Lf* 
Lp t 
Ll* 
L l . 
Lr* 
Lr, 
Lb* 
Lb* 

2/s 



2\Vu ~ V 2l , 7 r—E i R 
p + ie 



1 



2[U U - U 2l , 



1 



a + ie 



-El 



1 



2[U 1% -U 2i ,- 



1 



-El 



13 -ie 
1 

a — ie 
1 

~ie' 
1 



-El 



-El 



-El 



a + ie 

2[V U -V 2i ,-1—E% 
p + ie 

2[U U - U 2i , —^—E\ 



a — ie 
1 



-El 



-El 



2[V U - V 2ll - 7r --E F 



= [Uu-U 2l , 



where 7^±k°( x ) = 1 Jo ±0 ° du °i x + U P2) and -^-0{x) = 
i J ±o ° duO(x + upi). It is easy to check the background- 
Feynman gauge condition (id^ + g[W F ,)L Ftl = (and 
similarly for three other quadrants of the space). 

The transverse part Ei agrees with the results of Sec. 
while the longitudinal part l|A9(l does not literally agree 
with l|32l) (sec the footnote after that equation). It should 
be emphasized that, unlike the calculations with trial 
configuration (|25|l . the Feynman diagrams in the back- 
ground of the ansatz l|Alfl are free from uncertainties like 
0(0). 

Let us rederive now the effective action (|39[) starting 
from the ansatz (|A1|) and the fields (|A8|) . Since the only 
non-zero component of the field strength for the ansatz 
(|A1|) is transverse (see Eq. (|A3|) . we have 



a + ie 
1 

P + ie' 
a + ie 



-El 



13 -ie 
1 

a — ie ~~ 



-Eh] (A9) 



+ l - [ d i zd i z'T l a {z)THz')(A at {z)A h ^z')) 



= - - I d i z(Fj° )a F^ a ' ik 



+ i j d 4 z'F^ a (z)((&A ak (z) i - k)A^(z'))T b (z')} 
= -- A Jd A zF^ )a F^ k (A10) 



where F^ a,lk is a field strength in the first order in 
[U, V] . Using the fields l|A8|l we obtain 



(F L) + (F R) + (F B) - (i <-> k) 



F% )a {z) = -2g9(F)(z\wl^-W F \ ab 0,E b k ) (All) 
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(where \0,Ei) = J d 2 z' ± \0,z' ± )E i (z' ± )) and therefore 



S, 



oft 



(A12) 



ig 2 f <fiz{6{F){[UlV* 



-i^k) a {wl^W F ) ah E b k 



! a 2 
ft 
-a 2 
ft 



A typical integral in the above equation has the form 



* l<Pz x Mj(z x ){z±_\% W x Wx) (A13) 

27rJ a J p\ 



In the LLA, the integral J^°°^ is replaced by More 
accurately, one should remember that the slopes of Wil- 
son lines are e\ = Pi + e~ m p 2 and e 2 = P2 + e m pi as 
shown in Eq. (fHfl . In this case, 6(z*)9(z,) in the inte- 
grand of Eq. (|A13|) will be replaced by 6{z^+e r ' 2 z.)9{z,+ 
e~ m z*) so one obtains 



d'zSA 9(z*)9(z.)f(z x )(z\^-\0,z' ± )g(z' ± ) 



- d 2 z x d 2 z'J{z x )g{z' x ) 



P Z + IE 

dad(3d 2 p x i 
16tH 



(p,z-z')j 



a(5s - p\ + ie (a + e r ' 2 P-ie)((3 + e- i n a -ie) 
- i [ d 2 z x d 2 z' ± f(z x )g{z' x ) [°°-^ e t(p ' Z - Z ' }± 



a 2 s+p ± e , ' 2 a 2 s+p 2 



= --LAi, d 2 z ± d 2 z' ± f(z x )(z^\z' ± )g(z' ± ) (A14) 

where A?y = 771 — rj 2 . Performing the integrations over 
z* , z. in Eq. i|A12(> we get 



(A15) 



At? y ^{([CT^*] - i «-> k) a (wl^W F ) 
- ([Ui,Vfl-i~k) a (W&W R ) ab 



d 
ft 



([i^X]-^fc) a (^^-^) Qb 
([ui v 2 k ] i <- kYiw^WBY^E, 



= - ia s Ari J d 2 z x E?E ai 
which coincides with Eq. (|39|l . 




FIG. 11: Trace of four Wilson lines. 



Finally, let us demonstrate that the "diamond" trace 
of four (non-differentiated) Wilson lines is trivial (this is 

related to the fact that the field strength vanishes in 

the leading order, see Eqs. i|A8(l and (jA9|0 . To regularize 
the corresponding expressions, we consider the "original" 
tilted Wilson loop shown in Fig. ^2 for the finite size L. 
We need to prove that 

lim tr{[— Lei + Le 2 + x x ,Lei + Le 2 + x x ] (A16) 

x [Le 2 + Lei + x x , -Le 2 + Lei + x x ] 
x [Lei - Le 2 + x x , -Lei ~ Le 2 + x x ] 
x [—Le 2 - Le 2 + x x , Le 2 - Le 2 + x x }} = 1 

in the leading nontrivial order in [U, V] . 

Consider the case Ui <C 1, Vj ~ 1 (the opposite case 
Vi -C 1, Ui ~ 1 is similar). It is easy to see from Eq. 
(1A"9|) that [Le 2 ± Lei + x x ,-Le 2 ± Lei + x x ] ~ [U, V} 2 
so we are left with 

lim tr{[— Lei + Le 2 + x x , Lei + Le 2 + x x ] 

L — >oo 

x [Lei - Le 2 + x x , -Lei - Le 2 + x x }} (A17) 

At this point, we can take the limit L — * 00 in the ei 
direction. We obtain: 

[x x + Le 2 , x x + Le 2 + ooei] (A18) 
x [x x — Le 2 + ooex,x x — Le 2 ] — 1 



dadft 



sinaL 



x (x x \U} 



(a + e r ' 2 /3 - ie) (/3 + e^'' 1 a - ie) 
1 



a(3s — p~i_ + i 



-U 1 \ ab \[U li -U 2U El-El] b ) 



2 -r d a^(x x \uU i-^- 

irj a 1 \e-^a 2 s + p\ 



1 



e -»?2Q,2 s + p 2 ± 



Ui\ ab \[Uii-U 2i ,El-E^) 



We see now that in the limit L — > 00 the r.h.s. of Eq. 
(|A19(1 vanishes so the l.h.s. is at best ~ [U, V] 2 . Simi- 
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larly, 



and therefore the trace (|A17|I . which is product of l.h.s. 
of Eq. (|A18|) and Eq. I|A19|I , is equal to 1 in the leading 



lim [x±_ + Le2, x± + Le2 — ooe{\ 



(A19) order. 



x — Le2 — ooei, x± — Le^] = 1 
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